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INTRODUCTION 


We have completed our study of the traceability of the calibrated COsmic Background 
Explorer/Diffuse InfraRed Background Experiment (COBE/DIRBE) Photometric Standard Value 
(PSV) catalog to the "accepted" values as given in Gezari’s "Catalog of Infrared Observations" 
(CIO). This study encompasses nearly 100 celestial sources over a wide range of intensity and 
temperature. The results provide an independent statistical verification of the DIRBE point 
source calibration. 

The official COBE/DIRBE calibration is based on a single celestial source for each band. The 
current calibration, called "Mark 3", provides a transformation from linearized DIRBE digital 
numbers (DNs) to flux as described in the following table: 

DIRBE Mark 3 Calibration -11 July 1995 


Ban 

d# 

Source 

(w/m 2 -Sr^j 

Jy_ 

D N 

MJy/Sr 

DN 

ES^SRS 

uncertainty 

% 

MJy/Sr 

DN 

1 

Sirius 

2.46468e-007 

51.035 

0.426 

1.0224e-006 

3 

2 

Sirius 

8.86400e-008 

55.522 

0.391 

5.3318e-007 

3 

3 

Sirius 

7.56593e-008 


0.338 

2.897 le-007 

3 

4 

Sirius 

2.45934e-008 

43.597 

0.298 

1.8245e-007 

3 

5 

Sirius 

1.64531e-007 

161.108 

1.129 

2.8225e-007 

4 

6 

NGC7027 

6.15779e-008 

212.139 

1.457 

1. 7484c -007 

15 

7 

Uranus 

1.17295e-007 

722.736 

4.780 

2.3900e-007 

11 

1 8 

Uranus 

3.91026e-008 

556.605 

3.906 

1.1718e-007 

13 

9 


4.7982 le-007 

10539.851 

76.100 

1.6307e-006 

11 

10 


2.09268e-007 

5358.150 

40.500 

5.0625e-007 
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This study reports an independent statistical check on this calibration using an ensemble of 
celestial sources for each DIRBE band. 


2 DIRBE INSTRUMENT DESCRIPTIONS 

DIRBE is a 10 band, single element, superfluid He cooled, off-axis Gregorian telescope. 
Following the Mark3 conventions, the DIRBE wavebands are defined as: 


DIRBE Band 

Wavelength 

(microns) 

Bandwidth 

(microns) 

Beam Size 
(steradian) 

1 

1.25 

0.30 

1.98E-04 

2 

2.20 

0.37 

1.42E-04 

3 

3.50 

0.91 

1.29E-04 

4 


0.66 

1.46E-04 

5 

12.00 

7.00 

1.43E-04 


1 

























































































6 

25.00 

8.80 

1.46E-04 

7 

60.00 

29.45 

1.5 IE-04 

8 

100.00 

33.37 

1.43E-04 

9 

140.00 

41.19 

1.39E-04 

10 

240.00 

99.21 

1.32E-04 j 


3 METHODOLOGY 

We performed a blind test of the new Mark 3 calibrated DIRBE point source PSV catalog 
obtained from the COBE web site: 

http://bolero.gsfc.nasa.gov/astro/cobe_home.html 

The PSV catalog contains DIRBE-constructed PSVs in janskys for 92 BEX_COL (DIRBE 
External Calibrator Object List) sources in all bands for which the source was an approved 
DIRBE calibrator, and for which an average could be constructed. The catalog also contains the 
RMS, the backgrounds, and the number of observations. The PSV catalog contains sources 
useful for calibration from 1-200 microns, covering non-stellar objects and stellar spectral types 
from A to M over a wide dynamic range. The DIRBE PSV values are produced by averaging 
over typically 500 single observations of each source, and have been background-subtracted 
using spatial differencing. The typical RMS for the resulting averages are about 10 3 of the 
source flux. 

Our objective was to create a "predicted" PSV catalog from the CIO. We used the 1994 version 
of the Catalog of Infrared Observations (CIO) from the NSSDC (National Space Science Data 
Center). NASA Reference Publication 1 196, 3rd Edition, 1994, Gezari, Schmitz, Mead. The 
1994 version contains more years of data (both earlier and later) from the literature than does the 
1 987 ADC CD-ROM version. The CIO is an un-edited compilation of values from the literature, 
with native units spanning many systems from magnitudes, Janskys, W/cm 2 "inband", W/cm 2 
"XTk", radiance, irradiance, spectral, inband, and so on. 

We first filtered the CIO into a smaller catalog containing all observations of the 92 BEX_COL 
objects in bands that overlapped DIRBE bands. We then derived photometric transformations 
from the various CIO flux units into both janskys (the units of the DIRBE catalog) and W/(cm 2 - 
sr-pm). 


Because the CIO only contains effective wavelengths but not bandwidths, we were able to 
transform the native CIO observations to the proper effective wavelength, but were unable to 
correct for bandwidth for the CIO inband observations, as that would require manually 
researching the original CIO source literature which is beyond the scope of this study. No color 
corrections were attempted with the CIO data, as it would again require manually reviewing the 
CIO source literature. 
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We then created the "CIO predicted" PSV catalog by robust averaging all the ClO-observations 
of each BEX_COL source in each band, keeping the rms and number of observations. Extreme 
outliers in this distribution were discarded. 

We next did a catalog search for the spectral types of all the stellar BEX_COL sources. We were 
able to find spectral types for over 3/4 of the sources. We assigned an effective blackbody 
temperature to each star based on its spectral type. Finally we applied a DIRBE-band color 
correction to the calibrated DIRBE fluxes. Sources, particularly the non-stellar sources, that we 
could not get a spectral type or temperature for were not color corrected. 

As will be discussed more fully later, we used the Sharpless catalog to get the spatial size of 
several of the longer wavelength sources. Given the source sizes in conjunction with the DIRBE 
and CIO reported beam sized, we were able to do a primitive spatial extent correction on the 
appropriate non-stellar data. 

Finally, we compared the CIO-predicted and observed PS Vs, calculating RMS delta_flux, and 
flagging all those observations which did not fit into each others RMS (root-mean-square) 
uncertainties. We did this in both janskys and W/(cm 2 -sr-pm). We plotted a flux vs. flux 
correlation plot for all the sources in each DIRBE band. We calculated two statistical metrics: 
The correlation coefficients: 

k Zhh 


and the band-averaged percent differences: 

d = 



\ 

) 


Where the x* and y, are the DIRBE and CIO measurement sets respectively. There were typically 
~500 measurements in the DIRBE ensemble for each source for each band. The CIO in contrast 
typically contained only -10 measurements in an equivalent ensemble. 


4 RESULTS 

The enclosed plots and table in Appendices A & B respectively give the CIO-predicted and 
observed flux in W/(cm 2 -sr-|im) (abbreviated as "wcsu"), RMSs, number of observations, and 
DERBE-derived backgrounds for all CIO observations of DIRBE BEX_COL sources for all 
bands in which the source is an approved calibrator. Where sources or bands are missing it is 
because either the CIO did not contain relevant observations, or DIRBE was unable to construct a 
PSV. The final four columns contain the absolute and percent differences between the predicted 
and observed RMS PS Vs, as well as a flag for those RMS observations which did not fit into 


3 



each others RMS uncertainties (0:agree, l:disagree). The last column gives the assigned source 
spectral type where available. 

There was no overlapping data between the CIO and DIRBE for band 10. 

The figures in Appendix A are correlation plots between DIRBE and CIO flux. If DIRBE and the 
various CIO observations were all in perfect agreement, all the points would fall on the 1 : 1 line 
indicated. The CIO errorbars (horizontal) are from the RMS of the various CIO measurements of 
that star in that band. The DIRBE errorbars (vertical) are the DIRBE RMS error. The DIRBE 
errors are typically too small to show up on the plots, though they are tabulated below. 

The "band-average-percent difference" d is a measure of how well the DIRBE points agree with 
the CIO points on average, over each DIRBE band. 

The correlation coefficient k is a measure of how well the data fits the 1:1 line. 


5 CONCLUSIONS 


5A GENERAL CONCLUSIONS 

DIRBE was able to construct internal PSVs that agree with the literature, especially in the shorter 
wavelengths (bands 1-6). Longer wavelengths (bands 7-9) are not as good. Remember that we 
did not typically color correct the non-stellar fluxes. There are only 3 overlapping sources in 
each of the bands 7-9. and they are all extended sources. This leads to possible systematic errors 
that will be discussed later. 

The agreement between the predicted and observed PSV over many bands and many sources is 
encouraging. Almost all the disagreement (deltas of more ~5 %) we attribute to beam size 
effects due to bright backgrounds, or the source having extended IR structure. Some of the 
remaining anomalies may be related to unknown or uncorrected bandwidths. Many of the CIO 
observations are ground-based, so errors in the original investigators’ atmospheric corrections 
may also contribute to the differences. 

This current work represents a significant improvement over the preliminary version that we 
reported last year. The previous version used the older CIO database, the older "mark 2" DIRBE 
calibration, and performed no color or beam size corrections. With these improvements, the 
typical percent difference between the cio-derived predicted PSVs and the actual DIRBE PSVs 
has been significantly reduced. 
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5B SHORT WAVELENGTHS 


As expected, in DIRBE bands 1-4 (1-5 microns), we have many clean stellar point sources, 
typically with good temperature/spectral type corrections. The correlations and percent 
differences are excellent. 


5C MEDIUM WAVELENGTHS 

In DIRBE bands 5 and 6 (12 and 25 microns) we begin to have several non-stellar sources in the 
comparison. Because the zodiacal light is integrated over a large field of view, the point sources 
are typically 10-100 times fainter than the average background. The pipeline processing 
algorithm which generated the PSV catalog subtracted an average background over the 10 
months of observations. Because the Zodiacal Light is a strongly time dependent background, 
this leads to a residual systematic background error. This average background can be corrected 
in principle by re-creating a point-source average value from the real data, but this time 
subtracting the zodiacal signal for the appropriate position and date prior to averaging, but this 
effort is beyond the scope of this study. 

There are numerically fewer comparison sources in the statistical samples for the medium 
wavelength regime. Even with these limitations, the correlations and percent errors are still 
good. 


5D LONG WAVELENGTHS 

In the longer 60-140 micron wavelength bands (DIRBE bands 7-9) There are many possible 
systematic and random errors. One of the biggest limitations is the very small number of sources 
available: only 3 sources in each band. This makes the statistics inherently poor. In addition, 
every source in these bands is an extended cool object, typically much larger than the CIO beam 
sizes. The Sharpless catalog reports the angular sizes to be: 


Source 

Name 

Angular Size 
(arcmin) 

# Assoc 
Stars 

Typical CIO_beam 
Size (arcsec) 

S140 

30 

1 

120 

SI 84 

40 

1 

30 

S222 

6 

0 

37 

S235 

10 

0 

180 

S254-257 

~15 

1 

50 


Remember that the DIRBE beam size was about 1.5e-4 sr, or about 0.7 deg. 

In addition the Sharpless catalog indicates that all these sources are morphologically 
"elliptical/filamentary", i.e. exhibit strong spatial structure. We did only a rudimentary area-ratio 
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correction to the CIO flux, essentially calculating what the CIO measurement should have been if 
it spanned the entire spatial extent of the structure, and assuming that the source was spatially 
uniform. This is expected to be a gross overcorrection, so we assigned the resulting CIO 
uncertainties to be 100%. 

The band 7-9 plots show two sets of points: the beam-uncorrected data are the standard boxes 
with solid CIO errorbars. The beam-corrected data are the * symbols with the dashed CIO 
errorbars. Before the angular size corrections, the DIRBE data is systematically higher than the 
CIO data, but still roughly parallel to the ideal 1:1 correlation line. As expected, because the CIO 
beams are typically much smaller than the sources, while DIRBE typically integrates the entire 
spatial extent. 

Systematic effects are also expected due to the lack of color corrections in either these non-stellar 
objects or the backgrounds, though the spatial corrections dominate. 

It should be noted that the detectors in these long wavelength bands also exhibited extreme non- 
linear behavior, with uncorrected residuals typically near 10-20 percent. 
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Appendix A: 

Presentation given at session 51.03 of AAS 189, 14 January 1997, Toronto, Ontario, Canada. 
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Abstract 




We report the results of a comparative study 
of the COsmic Background Explorer /Diffuse 
InfraRed Background Experiment (COBE/DIRBE) 
photometric calibration over ~100 selected stellar and 
non-stellar calibration objects across a wide range of 
the DIRBE instrument dynamic range, wavelength 
coverage, and source temperature. A statistical 
comparison of the DIRBE-reported flux to the 
accepted values from the literature (as summarized 
in the CIO) provides an independent verification of 
the DIRBE point source calibration. 


9701A002 

aas197sb.ppt 
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The DIRBE Instrument 



CSP 

Center for Science and Phcrhxnenob%y 


• 10 band Radiometer 


Band 

Wavelength 

(microns) 

Bandwidth 

(microns) 

Beam Size 
(steradian) 

1 

1.25 

0.30 

1.98E-04 

2 

2.20 

0.37 

1.42E-04 

3 

3.50 

0.91 

1.29E-04 

4 

4.90 

0.66 

1.46E-04 

5 

12.00 

7.00 

1.43E-04 

6 

25.00 

8.80 

1.46E-04 

7 

60.00 

29.45 

1.51E-04 

8 

100.00 

33.37 

1.43E-04 

9 

140.00 

41.19 

1.39E-04 

10 

240.00 

99.21 

1.32E-04 


• Single element He cooled detectors 

• Off-axis Gregorian telescope 


9701A003 
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Methodology: 
Overview 

• We compare the DIRBE observed flux to the accepted literature: 

- Wide wavelength range (1-200 microns) 

- Wide temperature range (A-M spectral type stars...) 

- Multiple sources (~100) 

- Wide dynamic range 

• DIRBE: 

- The 92 source DIRBE Photometric Standard Values (PSV) obtained 
from the COBE web site: 

» http://bolero.gsfc.nasa.gov/astro/cobe/cobe_home.html 
» Typically 

- ~ 500 observations of a given calibration object 

- Background subtracted and averaged 

- RMS ~ 10' 3 of source flux 

• Accepted Values: 

- The Catalog of Infrared Observations (CIO), NASA Reference 
Publication 1196, 3rd Edition, 1994. Gezari, Schmitz, Mead. 




9701A004 
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Methodology: 
DIRBE Processing 



• Transform the DIRBE catalog entries from Janskys to 
W/ (cm 2 -sr-jim). 

• Look up spectral type and/or temperatures for most DIRBE 
sources. 

• Apply DIRBE instrument color corrections to observations. 


9701A005 

aasi97sb.ppt 
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Methodology: 
CIO Processing 



• Filter the CIO electronic catalog for all entries containing DIRBE 
PSV sources that overlap the DIRBE wavelength bands. 

• Transform the selected CIO entries to W/(cm 2 -sr-|im). 

• Calculate weighted robust average & RMS of all CIO values: 

- For each DIRBE source 

- For each DIRBE wavelength 

• Create a "predicted" DIRBE PSV catalog subset from resulting 
CIO database. 

• No color corrections done with CIO data. 

- Not possible without going to original source literature 

- Little expected gain for large effort 


9701A006 
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Methodology: 

Comparison 



• Correlation plots by DIRBE band for overlapping data. 

• Error bars in both DIRBE & CIO axes from the RMS of 
multiple observations. 

• Calculate correlation coefficients: 

k=sum(x i *y i ) / sqrt(sum(x i 2 )*sum(y i 2 )) 

• Calculate band-average percent errors: 

d^sumCabsCx-yj))/ (sum(x i +y } )/2) 


9701A007 

aas197sb.ppt 
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Cautions 




• Possible systematic errors for longer wavelengths: 

- Temperature/ color correction errors: 

» Background 
» Source 

- Beam size effects: 

» Extended sources 
» Background effects 

• CIO is a "metacatalog" of mixed quality. 


9701 A008 
aas197sb.ppt 
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Discussion: 
Shorter Wavelengths 



• DIRBE bands 1-4 (1.25 -4.9 pm): 

- Stellar sources: 

» Point-like sources 
» Good temperature/ spectral type data 

- Many sources per waveband: good statistics. 

- Excellent correlation & good calibration. 


9701 A010 
aasl97sb.ppt 
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Discussion: 
Medium Wavelengths 



• DIRBE bands 5 & 6 (12 - 25 jum): 

- Some non-stellar sources 

- Sources typically 10-100 times fainter than background 

- Not as many common DIRBE /CIO sources 

- Good correlations & reasonable calibration 


9701 A011 
aas197sb.ppt 
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Discussion: 
Longer Wavelengths 



• DIRBE bands 7-9 (60 - 140 pm): 

- Few overlapping sources: poor statistics 

- Roughly parallel with ideal 1:1 line 

- DIRBE systematically higher than CIO values 

- Systematic errors: 

» Temperature/ color corrections 
» Beam size /source size 

• DIRBE band 10 (240 jum): 

- No overlapping data 


9701A012 

aas197sb.ppt 
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Band 7, 8, 


9 Systematics 



• 6 sources: 5 Sharpless HII regions & M82. 

• Estimate beam size correction factors & errors: 

- 42 min DIRBE beam vs typically -30-120 sec CIO beams 

- Sharpless sources typically: 

» 10-50m angular size 
» Elliptical, filamentary shape 

- Implies DIRBE values up to -50 times greater than CIO 

- The correction factors are upper bounds 

- Apply Corrections to CIO data: 

» Plotted as * symbols 

» Corrected CIO uncertainty taken as 100% (dashed lines) 


9701A013 

aas197sb.ppt 
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Appendix B: 


The Following table is the result of the IDL analysis routines used to process the data in this 
study. The columns contain: 


bex# 

bex_name 

band 

CIO_wcsu 

CIO_rms 

CIO_#obs 

psv_wcsu 

psv_rms 

psv_bkg 

psv_#obs 

delta 

ffac 

flag 

spctyp 


DIRBE BEX_COL (External Calibrator Object List) number 
DIRBE BEX_COL source name 
DIRBE band 

CIO average flux in w/(cm 2 -sr-pm) 

CIO rms uncertainty 

number of observations in CIO ensemble for this source in this band 
DIRBE average flux in w/(cm 2 -sr- pm) 

DIRBE rms uncertainty 

DIRBE constructed background flux 

number of observations in DIRBE ensemble for this source in this band 
the difference between the CIO and DIRBE flux in w/(cm 2 -sr- pm) 
the ratio between the CIO and DIRBE flux 
a flag to indicate top-level agreement: 0: agree, 1: disagree 
source spectral type if known 
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Detailed Table of DIRBE/CIO Source Comparison Results: 

Comparison of DIRBE calibrated PSC vs CIO average values: run date: 29-Dec- 1996 20:42:29.00 
Predicted DIRBE JYs (avg of CIO BEX_COL obs w/ mark3 cal) GRC/burdick 19-Dec- 1996 


CIO 


bex# bex_name 

band 

CIO_wcsu 

CIO_rms 

_#obs 

psv_wcsu 

psv_rms 

psvbkg 

psv_#obs 

delta 

frac 

flag 

spctyp 

27 

ALP-CAS 

1 

1.76e-009 

0.00e+000 

2 

1.812e-009 

1 ,476e-01 2 

1.408e-010 

674 

2.230e-01 1 

0.012 

1 

K0 

27 

ALP-CAS 

2 

3.98e-010 

3.74e-011 

3 

3.476e-010 

2.555e-013 

3.620e-01 1 

613 

5.491e-01 1 

0.160 

1 

K0 

27 

ALP-CAS 

3 

9.29e-01 1 

0.00e+000 

2 

7.21 le-01 1 

5.818e-014 

8.900e-012 

654 

2.059e-01 1 

0.285 

1 

K0 

27 

ALP-CAS 

4 

1.04e-01 1 

0.00e+000 

1 

1.499e-01 1 

1.630e-014 

7.127e-012 

716 

4.823e-012 

0.316 

1 

K0 

115 

S184 

7 

1.49e-013 

5.70e-014 

2 

2.012e-012 

9.121 e-0 1 5 

1.075e-012 

342 

2.753e-012 

0.949 

1 

M0 

115 

S184 

8 

1.59e-013 

0.00e+000 

1 

1.892e-012 

9.233e-015 

6.476e-013 

300 

2.104e-012 

0.930 

1 

MO 

71 

BET-AND 

1 

5.36e-009 

1.54e-009 

14 

5.996e-009 

3.992e-012 

6.808e-01 1 

364 

5.602e-010 

0.095 

0 

MO 

71 

BET-AND 

2 

1.57e-009 

3.66e-010 

23 

1.553e-009 

5.295e-013 

1.563e-01 1 

366 

3.372e-01 1 

0.022 

0 

MO 

71 

BET-AND 

3 

3.88e-010 

2.84e-01 1 

15 

3.353e-010 

1.486e-013 

3.998e-012 

350 

5.23 le-01 1 

0.156 

1 

MO 

71 

BET-AND 

4 

7.44c-011 

7.95e-012 

14 

6.458e-01 1 

2.55 le-01 4 

6.424e-0I2 

463 

8.745e-012 

0.133 

1 

MO 

71 

BET-AND 

5 

3.04e-012 

3.01e-013 

38 

2.588e-012 

1.229e-014 

3.693e-01 1 

345 

4.372e-014 

0.015 

0 

MO 

88 

GAM-AND 

1 

2.61e-009 

1.37e-010 

4 

2.776e-009 

1.460e-012 

8.250e-01 1 

479 

1.239e-010 

0.045 

0 

K3 

88 

GAM-AND 

2 

6.38e-010 

2.76e-01 1 

5 

5.963e-010 

2.141e-013 

1.912e-01 1 

481 

4.815e-011 

0.082 

1 

K3 

88 

GAM-AND 

3 

1.34e-010 

1.01e-011 

6 

1.276e-010 

5.896e-014 

4.678e-012 

488 

6.440e-012 

0.050 

0 

K3 

88 

GAM-AND 

4 

2.52e-01 1 

0.00e+000 

1 

2.505e-011 

1.666e-014 

6.263e-012 

573 

2.326e-013 

0.009 

1 

K3 

113 

ALP-CET 

1 

4.61e-009 

3.03e-010 

11 

4.978e-009 

3.425e-012 

6.707e-01 1 

346 

3.117e-010 

0.063 

1 

M2 

113 

ALP-CET 

2 

1.39e-009 

7.28e-01 1 

13 

1 .3 1 8e-009 

4.535e-013 

1.458e-01 1 

340 

8.514e-01 1 

0.065 

1 

M2 

113 

ALP-CET 

3 

3.19e-010 

2.53e-01 1 

17 

2.885e-010 

1.596e-013 

4.139e-012 

338 

3.033e-011 

0.105 

1 

M2 

113 

ALP-CET 

4 

6.24e-01 1 

6.21e-012 

8 

5.469e-01 1 

2.516e-014 

7.022e-012 

432 

6.874e-012 

0.124 

1 

M2 

93 

RHO-PER 

5 

3.41e-012 

6.65e-013 

9 

3.142e-012 

1.162e-014 

3.998e-01 1 

419 

2.08 le-01 3 

0.057 

0 

M4 

120 

S222 

3 

4.15e-01 1 

7.15e-012 

5 

4.392e-01 1 

8.019e-014 

7.730e-012 

255 

2.392e-012 

0.054 

0 


120 

S222 

4 

2.59e-01 1 

3.04e-012 

4 

3.000e-01 1 

4.678e-014 

7.753e-012 

355 

4.06 le-01 2 

0.135 

1 


120 

S222 

8 

4.60e-013 

9.09e-013 

12 

3.708e-012 

1.015e-014 

6.880e-013 

163 

3.247e-012 

0.876 

1 


112 

ALP-TAU 

1 

1.49e-008 

1.03e-009 
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